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OO \ In this article we study from a non-perturbative point of view the entanglement of 

^ ^ two directed polymers subjected to repulsive interactions given by a Dirac <5— function 

^ ■ potential. An exact formula of the so-called second moment of the winding angle is 

. derived. This result is used to provide a thorough analysis of entanglement phe- 
nomena in the classical system of two polymers subjected to repulsive interactions 

C ' and related problems. No approximation is made in treating the constraint on the 

^ ■ winding angle and the repulsive forces. In particular, we investigate how repulsive 

ff^ \ forces influence the entanglement degree of the two-polymer system. In the limit of 
> 

'sj" ' ideal polymers, in which the interactions are switched off, we show that our results 

^ ■ are in agreement with those of previous works. 

o 



I. INTRODUCTION 



pH ■ The statistical mechanics of two polymers with constraints on their winding angle has 

O '. 

O . been extensively studied in order to understand the behavior of physical polymer systems, 



like for instance biological macromolecules of DNA L, 



or liquid crystals composed of stacks 



j_j ■ of disk-shaped molecules see Refs. 

detailed review on the subject, together with interesting proposals of how to include in the 
treatment of topologically entangled polymer link invariants which are more sophisticated 
than the winding number, can be found in Up to now, however, despite many efforts. 
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mainly ideal polymer chains or loops winding around each others have been considered, 
while the repulsive interactions between the monomers have been treated approximatively 
or exploiting in a clever way scaling arguments integrated by numerical simulations, as for 
instance in 

Here we concentrate ourselves on the case of two directed polymers interacting via a 
repulsive Dirac 5— function potential 0, [sijl- We are particularly interested in the average 
degree of entanglement of the system, which we wish to estimate by computing the square 
average winding angle of the two polymers. This quantity is also called second moment 
of the winding angle or simply second moment and is a special example of the topological 
moments first introduced in Ref. To achieve our goals, we develop an approach, which 
combines quantum mechanical and field theoretical techniques. With respect to previous 
works, we are able to obtain exact results even if the repulsive interactions are not switched 
off. 

In principle, the average of any observable like the squared winding angle can be derived 
once the partition function of the system is known, but in our case it turns out that the 
partition function is simply too complicated to obtain any analytical result. This happens 
essentially because the full 5— function potential is not a central potential, since it mixes 
both radial and angular variables. For this reason, the usual procedure of going to polar 
coordinates and then solving the differential equation satisfied by the partition function 
of the entangled polymers with the method of separation of variables (l^, does no longer 
produce simple formulas as in the situations in which only central forces are present. 

To avoid these difficulties, one possibility is to approximate the 5— function potential 
with some radial potential, like for instance the hard core potential of Ref. However, 
here we shall adopt a different strategy, based on field theories, which does not require any 
approximation. This strategy has been RQ,ee also Kef. Q for more 

details) to cope with ideal closed polymers whose trajectories are concatenated. Also such 
systems are characterized by a non-central potential, which comes out as a consequence of 
the topological constraints imposed on the trajectories. In the field theoretical formulation 
of Refs. 0, [3, [l^ the computation of the second moment is reduced to the problem of 
computing some correlation functions of a field theory. A bonus is provided by the fact that 
this computation requires just a finite number of Feynman diagrams to be evaluated. In 
the present case, due to the presence of the 5— function potential, the field theory which we 



obtain is no longer free as that of Refs. 



E, 12, 3 



Nevertheless, we will see that the theory 



is still linear and thus it can be exactly solved once its propagator is known. Luckily, this 
propagator may be computed exactly using powerful non-perturbative techniques developed 
in the context of quantum mechanics to deal with Hamiltonian containing 5— function po- 
3, 24,0, 2£, 22,0,0,0,0,0- 

Basically, starting from the Green 



tentials, see Refs. 



function of a particle whose dynamics is governed by a given Hamiltonian Tio, these tech- 
niques provide an algorithm to construct the Green function of a particle corresponding to 
a perturbed Hamiltonian Ti. = Tio + Vs, where Vs is the 5— function potential. One advantage 
of these methods is that there is a long list of potentials for which the Green functions of the 
unperturbed Hamiltonians Hq are known. In this way, it is easy to generalize our treatment 
including new interactions, which could be relevant in polymer physics, like for instance the 
Coulomb interaction. The price to be paid is that the quantum-mechanical algorithm works 
when the Green functions are expressed as functions of the energy instead of the time. In 
the polymer analogy, assuming that the ends of the polymers are attached to two planar 
surfaces perpendicular to the 2;— axis and located at the positions z = and z = L, the role 
of time is played by the distance L, while the energy corresponds to the chemical potential 
conjugated L. To recover the original dependence on L, one needs to calculate an inverse 
Laplace transform of the field propagator with respect to the energy. In general, this is not 
a simple task. 

Once the propagator of the linear field theory is known, the correlation functions which 
enter in the expression of the second moment may be calculated contracting the fields in 
all possible ways using the Wick prescription. At the end, we get in this way an exact 
formula of the second moment as a function of the energy, which, we remember, has here 
the meaning of the chemical potential conjugated to the distance L. In the L space, due 
to the problems of computing the inverse Laplace transform of the propagator mentioned 
above, only an approximated expression of the second moment will be given in the limit of 
large values of L and assuming that the strength of the 5— function potential is weak enough 
to allow a perturbative approach. 

Our results allow both a qualitative and quantitative understanding of the way in which 
the repulsive interactions affect the entanglement of two directed polymers. The corrections 
introduced by these interactions in the expression of the second moment of ideal polymers 
have been studied in some interesting limits. First of all, it has been examined the limit 
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of long polymer trajectories, in which we show that repulsive interaction become particu- 
larly relevant. Moreover, we have investigated also the perturbative regime and the strong 
coupling limit, which is important to recover the excluded volume interactions. While it is 
not a problem to take the strong coupling limit within our exact treatment of the repulsive 
interactions, it turns out that, in this case, the expression of the second moment is partic- 
ularly complicated from the analytical point of view. For this reason, in the Conclusions 
we will discuss the application of a powerful perturbative method to study field theories at 
strong coupling due to Kleinert 2^,13, 0|- Finally, the consistency of our results with the 
previous ones has been checked by studying the limit of ideal polymers. 

The material presented in this paper is divided as follows. In the next Section, the 
problem of computing the second moment of the winding angle of two directed polymers 
interacting via a 5— function potential is briefiy illustrated using the path integral approach. 
A constraint on the winding angle is imposed by coupling the trajectories of the polymers 
with a suitable external magnetic field, following the strategy of previous works like for 
instance j^, lU, Q, In Section ITTTl the second moment is expressed in the form of a 

finite sum of amplitudes of a linear field theory. These amplitudes may be computed once 
the propagator of the theory is constructed. In our case, the propagator coincides with 
the Green function of a particle diffusing in a 5— function potential. The derivation of this 
Green function in the energy representation using non-perturbative techniques developed in 

the context of quantum mechanics y , Q, Q, 0, H Q, Q, 0, Q is the subject of 



Section Hvl The (5— function potential is responsible of the appearance of singularities in the 
propagator at short distances, which have been regulated here with the introduction of a 
cut-off. This procedure is motivated by the fact that in polymer physics there is no point in 
considering distances which are smaller than the dimensions of a monomer. A comparison 
with the more rigorous method of renormalization is made, showing the consistency of 
the two procedures. The propagator derived in Section (llVjl has a particularly nice form, 
in which the contributions coming from the repulsive forces can be separated from the free 
part of the propagator, which is related to the random walk of ideal polymers. This splitting 
of the propagator is used in Section El to discuss qualitatively and qualitatively the effects 
of the (5— function interactions on the entanglement of the system. The results of Sections 
TTTl and IIVI provide in principle all the ingredients of the second moment. However, the 
amplitudes of the linear field theory derived in Section IIIII should still be evaluated. In 
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this task one encounters the typical problems occurring in the evaluation of the analytical 
expressions of Feynman diagrams. In the case of the second moment there are just tree 
diagrams, but still one has to perform complicated integrations over the spatial coordinates 
which are transverse to the 2;— axis. Even assuming that polymers are ideal, the analytical 
evaluation of these integrations requires drastic approximations, see for instance ^|. To 
avoid these difficulties, we average the second moment with respect to the positions of the 
endpoints of the two polymers. This averaged version of the second moment can be computed 
without any approximation in the energy representation. This is done in Section IVll The 
expression of the averaged second moment in the L— space is provided instead only at the 
first perturbative order in the strength of the repulsive potential and assuming additionally 
that the value of L is large. We give also an exact formula of the second moment without 
performing any averaging procedure as a function of L. This formula is however explicit only 
up to the calculation of the inverse Laplace transform of the propagator derived in Section 
IIVI In Section IVIII we consider the situation in which the polymers are not interacting in 
order to allow the comparison with previous results. Finally, the discussion of the obtained 
results and ideas for further developments are presented in the Conclusions. 



II. THE STATISTICAL MECHANICS OF TWO DIRECTED POLYMERS WITH 

CONSTRAINED WINDING ANGLE 



Our starting point is the action of two directed polymers A and B 

2 / , \ 2 



Ao = dz 
Jo 



dz 



dz 



(1) 



where V{va — ^b) is the potential: 



V{rA - tb) = -vo6{rA - vb) 



vo>0 



(2) 



The sign of Vq has been chosen in such a way that the interaction associated to the potential 
V{r) is repulsive. The parameters c and L determine the average length of the trajectories 
of the polymers. The ends of the polymers are supposed to be fixed on two surfaces perpen- 
dicular to the -z— axis and located at the heights z = and z = L. Both polymers have a 
preferred direction along the z direction. The vectors va^z) and rB^z)^ < z < measure 
the polymer displacement along the remaining two directions of the space. 
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The action of Eq. (QJ) resembles that of two quantum particles in the case of imaginary 
time z. To stress these analogies with quantum mechanics, the 2r— variable will be treated 
as a pseudo-time and renamed using from now on the letter t instead of z. 

In the system of the center of mass: 



va-^b R = (3) 



the action JT]) becomes: 



Ao= [ dt 
Jo 



(4) 



The motion of the center of mass, which is a free motion described by the coordinate R(t), 
will be ignored. 

We consider the partition function of the above two-polymer system with the addition of 
a constraint on the entanglement of the trajectories: 

Z^ = J Pre-/.^*[t(t)+v^M]5(^ - ^) (5) 
X is the so-called winding angle. Its expression is given by: 

X= rA(r(t))-dr(t) (6) 



where A(r) is a vector potential with components: 

In the above equation we have represented the vector r using cartesian coordinates 
i. e. r = (a;^,x^). Moreover, from now on, middle latin indices . . . = 1,2 will label the 
directions which are perpendicular to the t— axis. The definition of the partition function 
Zm is completed by the boundary conditions at t = and t = L: 

r(0) = ro r(L) = n (8) 

The quantity in Eq. © becomes a topological invariant if the polymer trajectories are closed. 
In the present case, in which the trajectories are open, x just counts the angle with which 
one polymer winds up around the other. Thus, the partition function Zm gives the formation 
probability of polymer paths winding up of an angle 

AO = 27rm (9) 
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Exploiting the Fourier representation of Dirac 5— functions 



/~f"CxD /V \ 
^Hmx) (10) 
-oo ZTT 



Eq. ^ can be rewritten as follows: 

^e-'Z, (11) 

where 

Zx = j Vre- (12) 

The Lagrangian £ is that of a particle immersed in the magnetic potential associated to the 
vector field (0): 

The Fourier transformed partition function Z\ is the grand canonical version of the original 
partition function Z^^ in which the number m is allowed to take any possible value. 

Zx coincides with the propagator Q\{L]Vi,Vq), which satisfies the following pseudo- 
Schrodinger equation: 

gxmr,,ro) = (14) 
H is the Hamiltonian of the system, computed starting from the Lagrangian (fT3|l : 

n = ^{V -tXAf + ViT) (15) 

Eq. is completed by the boundary condition at L = 0: 

e?A(0;ri,ro) = 5(ri-ro) (16) 

The average degree of entanglement of the two polymers can be estimated computing the 
topological moments of the winding angle (m^^)ri,ro) k = 0,1,2, .. . [2^. Once the partition 
function is known, the (m^^)ri,ro be expressed as follows: 

dm m^^Z^ dm m^^ f e-^g,(L; n, rp) 

The quantities (m^'^)ri,ro depend on the boundary conditions ro,ri and, of course, on the 
parameters c and L. For practical reasons, we will also consider the following averaged 
topological moments: 
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As Eq. (fT8|l shows, the average is performed with respect to the relative positions ro, ri of the 
endpoints. This is equivalent to an average over the positions of the endpoints rA(t),rBit) 
at the instants t = and t = L, because the coordinates of the center of mass have been 
factored out from the partition function and thus they do not play any role. The advantage of 
the averaged topological moments is that, a posteriori, it will be seen that their computation 
is easier than that of the topological moments given in Eq. (fT7|) . From the physical point of 
view, the averaged topological moments measure the entanglement of two polymers, whose 
ends at the instants t = and t = L are free to move. 

Here we will be interested only in the second moment {m'^)ri,ro the averaged second 

moment (m^), i. e. in the case A; = 1 of Eqs. (fT7j) and (f11?| . The second moment is in fact 
enough in order to estimate the formation probability of entanglement with a given winding 
angle and to determine how the winding angle grows with increasing polymer lengths. 

In the following it will be useful to work in the so-called energy representation, i. .e 
considering the Laplace transformed of the partition function Qx{L;ri,rQ) with respect to 
L: 



The new partition function Q\{E; ri, Vq) describes the probability of two entangled polymers 
of any length subjected to the condition that the relative positions of the polymer end at 
the initial and final instants to and ti are given by the vectors Vq and ri. With respect to the 
formulation in the L— space, however, the distance ti — to is no longer exactly equal to L, but 
is allowed to vary according to a distribution which is governed by the Boltzmann-like factor 
e^^. Thus, E plays the role of the chemical potential conjugated to the end-to-end distance of 
the polymer trajectories in the t— direction. It is worth to remember that, roughly speaking, 
small values of E correspond to large values of L, while large values of E correspond to 
small values of L. Starting from Eq. (fT^ and recalling the boundary conditions (fT6|) . it is 
easy to check that Qx{E]ri,rQ) satisfies the stationary pseudo-Schrodinger equation: 




(19) 



[E-n] ^?A(^;ri,ro) =<5(ri-ro) 



(20) 



where H is always the Hamiltonian of Eq. (fT5|l . 
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III. CALCULATION OF THE SECOND MOMENT USING THE FIELD 

THEORETICAL FORMULATION 

In this Section we wish to evaluate the expression of the second moment as a function 
of the energy E using a field theoretical formulation of the polymer partition function. The 
starting point is provided by the formula of the second moment in the L— space suitably 
rewritten in the following way: 

For consistency with Eq. (fTTjl . the numerator N{L; ri, Tq) and the denominator D{L; Vi, Tq) 
appearing in Eq. (f2T|l must be of the form: 

/+00 r+QO fj\ 

dmrn" — e^'"^6?A(L;ri,ro) (22) 
-oo J —CO ZTT 

and 

Dm n, ro) = / dm e'-^^g^^L- n, ro) (23) 

J —oo J —oo ZTT 

Using Eq. (fTHl) . it is now straightforward to compute the Laplace transform of N{L] ri, tq) 
and D{L; ri, tq): 

/ + 00 f + OO rj\ 

dmrn" —d^^g^{E-vr,v^) (24) 
-CxD J OO 2/K 

/+00 /'+00 A\ . 

dm\ —e'^^gx{E;r,,ro) (25) 
-oo ^oo ZTT 



Once the functions N{E;ri,ro) and D{E;ri,ro) are known, one can construct the ratio: 

which is nothing but the second moment of the winding angle expressed as a function of the 
chemical potential E. 

We remark that the Green function gx{E] ri,ro) is related to the Feynman propagator 
of the spin— I Aharonov-Bohm problem in quantum mechanics. In principle, this Green 
function can be computed exactly starting from Eq. (f2n|l but its final expression is 
too complicated for our purposes. Moreover, the method used in [3| to renormalize the 
singularities coming from the presence of the 5— function potential is valid only in a restricted 
region of the domain of A. This is incompatible with our requirements, because, to derive 
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the second moment, one has to integrate Gx{E; ri,ro) with respect to A over the whole real 
line. For this reason, we prefer here to use a field theoretical representation of this Green 
function. This is achieved by noting that Qx{E]ri,ro) coincides with the inverse matrix 
element of the operator E — Ti: 

g,{E-v,,To) = {v,\^^\vo) (27) 

and may be expressed in a functional integral form in terms of replica fields: 

GxiE; n, ro) = lim / I?vI/I?vl/*^,(ri)^*(ro)e~^('^*'*) (28) 

In the above equation \1/ are multiplets of replica fields: 

vl/ = (^i,...,^„) ^* = (^t, •••,<) (29) 

with action 



E-^^{V,-tXAf-vo6{x) 



^ (30) 



The symbol -k in Eq. ()28ll denotes summation over the replica index. For example t^t \1/ = 
I]ct=i V'ffV'cr- Below it will be used also the convention * = |\l>p. The details of the 
derivation of Eq. can be found in previous publications on the subject [2^ and will 
not be provided here. 

In order to proceed, it will be convenient to expand the action (f30|) in powers of A: 

^(^*, ^) = ^o(^*, ^) + A^i(^*, ^) + A2^2(**, ^) (31) 

where we have put: 

^o(^*,^) = Jd^x ^|V^|2 + (E-t;o5(x))|^|2 (32) 

5i(^*, xiJ) = l-J d^xA ■ i< (V^) - (V^*) ^ ^] (33) 

^2(^*, ^) = Ycj ^'2;A2|^|2 (34) 

At this point we come back to the computation of the quantities A^(i?;ri,ro) and 
D[E] ri, ro) appearing in the expression of the second moment. Exploiting the new form of 
the partition function given by Eqs. (f^HHHij) . together with the relation 

dm mV"^ = 2Ti{iY—^ 1/ = 0, 1, . . . (35) 
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and the fact that = 0, it is possible to rewrite Eqs. (pHl and (f25|l as follows |44l |: 

iV(E;ri,ro) = lim / PM/W^i(ri)^t(ro)[252(^*, ^) - (5i(vl/*, M/))2]e-^°(**'*) (36) 

D{E, n, ro) = lim / D^*I?^^i(ri)V'i*(ro)e-^«(**'*) (37) 

The right hand sides of Eqs. (p^ and (|H7jl represent vacuum expectation values of a field 
theory governed by the action S'o(\E'*, of Eq. In the formulation in terms of quantum 

operators we have: 

iV(E;ri,ro) = lim(0|^i(ri)^*(ro)252(^*, ^)|0)„-lim(0|V'i(ri)^i*(ro)(5i(^*, ^))2|0)„ (38) 

D(E;ri,ro) = lim(0|^i(ri)V^t(ro)|0)„ (39) 

The correlation functions have a subscript n to remember that, according to the replica 
method, they should be computed first assuming that the number of replicas n is an arbitrary 
positive integer and then taking the limit for n going to zero. 

The above correlators may be evaluated using standard field theoretical methods. One 
could be tempted to use a perturbative approach assuming that the value of vq appearing in 
the action 5*0 (^E^*, \E') of Eqs. and (jTrTjl is small, but this is not necessary. As a matter of 
fact, if it is true that S'o(\E'*, does not describe free fields because of the presence of the 
5— function potential, it is also true that it is just quadratic in the fields. As a consequence, 
one is allowed to define a propagator G{E; x, y) associated with this action. It is easy to 
check that G{E] x, y) satisfies the equation: 

1 



G(i?;x,y) = 5(x,y) (40) 



Using the above propagator, one can evaluate the amplitudes in Eqs. (j!^ and ex- 
actly by contracting the fields in all possible ways according to the Wick theorem. After 
straightforward calculations, one finds: 

lim(0|Vi(ri)V^i*(ro)|0)„ = G(E;ri,ro) (41) 

n— >u 

lim(0|V'i(ri)^r(ro)^2(^*,^)|0)„ = K{v,,r,) (42) 
lim(0|^i(ri)7/>*(ro)(5i(vl/*, ^)f\Q)n = h{r,, ro) + /2(ri, ro) + /3(ri, ro) + h{r,, ro) (43) 

n— >u 
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where 

Kin, ro) = ^ 1 d\A\^)G{E; n, x)G(E; x, ro) (44) 

/i(ri,ro) = |dW?/[A(x)G(E;x,ri)(V^G(E;y,x))A,(y)(V^yG(E;ro,y))](^^^ 

/2(ri,ro) = +^ |rf2xA[A,(x)(V5,G(E;x,ri))G(i5;;y,x)A,(y)(V^yG(E;ro,^ 

/3(ri,ro) = ciW|/[A,(x)G'(E;ro,x)(V^V^yG(E;x,y))A,(y)G(E;y,ri)] (47) 

/4(ri,ro) = |rfWy[Ai(x)G(E;ro,x)(V^G(E;x,y))A,(y)(V^yG'(E;y,ri))](48) 

From the physical point of view, the above equations may be interpreted in the following 
way. The fields \I/(x) and \l/*(x) contain operators which, inside each replica sector, create 
and annihilate segments of the two polymers, whose relative positions are given by the vector 
X. The two polymer system has been projected in the two-dimensional plane perpendicular 
to the t— axis. For this reason, there appear only the transverse coordinates x. The only 
remnant of the third dimension is the dependence on the energy E. The correlation functions 
(|i2|l and (|i3|) describe the fluctuations of the two polymers immersed in the 5— function 
potential and subjected to the interactions represented by the vector potential ((Zj). We 
recall that the origin of the latter interactions is the presence of the constraint on the 
winding angle in the partition function To evaluate the correlation functions (jl^ and 
(|^3j) . one needs to consider only a finite number of Feynman diagrams, corresponding to 
the relevant processes with which the two polymers interact together. The result, after the 
analytical evaluation of these diagrams, is provided by Eqs. (|i5l - li8|l . Let us note that in 
these equations the repulsive interactions due to the (5— function potential are hidden in the 
propagators (j(-E;x, y). The Feynman diagrams related to the amplitudes of Eqs. 
(jlH| are all topologically equivalent to the diagram of Fig. [H The amplitude of Eq. (jiill is 
related instead to the Feynman diagram of Fig. [2l The vectors ri and Tq denote the relative 
positions of the end points of the two polymers at the initial and final instants, as already 
mentioned. The integration variables x and y appearing in Eqs. (|iil - li8|l may be regarded 
as the vectors which give the relative positions of the trajectories of the two polymers at 
the instants in which they interact together via the external vector potential A of Eq. (0). 
There is no restriction on the domain of integration of x and y, so that the components of 
these relative vectors are allowed to take any value. This implies that the distance between 
the polymer segments when the interaction with A occurs can be arbitrarily large. 
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A 



A 



r 







FIG. 1: Feynman diagram corresponding to the amplitudes of Eqs. The two polymers 

A and B start at a distance |ro| from each other and interact twice with the the external field 
A. At the end, the relative position of the end points at the instant t = L is given by ri. The 
three- vertices appearing in the Figure are related to the interaction described by Eq. 



FIG. 2: Feynman diagram corresponding to the amplitude of Eq. H44|l . The two polymers A and 
B start at a distance |ro| from each other and interact with the the external field A. At the final 
instant t = L, the relative position of the end points is given by ri. The four- vertex appearing in 
this Figure is related to the interaction described by Eq. 

Now that the correlation functions which are present in the expressions of A^(£';ri,ro) 
and D(i?;ri,ro) given in Eqs. (|38|) and (j39|l have been evaluated, see Eqs. ()iil - li8|) . we may 
put everything together and give to the second moment of Eq. ()26p a more explicit form: 



In conclusion, the initial problem of computing the second moment of the winding angle 
(m^)ri,ro has been reduced to the evaluation of a finite number of integrals, which are given 
in Eqs. (|iil - li8|) . Of course, to make these integrals really explicit, we still need to derive the 
propagator G{E; x, y), which is so far the only missing ingredient. This will be done in the 
next Section. 



A 



A 



r 








(49) 
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IV. GREEN FUNCTIONS IN THE CASE OF HAMILTONIANS WITH A 

^-FUNCTION POTENTIAL 

Let Go{L; x, y) be the solution of the differential equation: 

(^^-^o^G'o(L;x,y) = (50) 

for a given Hamiltonian Hq. When L = 0, Go{L] x,y) satisfies the boundary condition: 

G'o(0;x,y) = 5(x-y) (51) 

In the case of a Hamiltonian H, obtained by adding to Hq a 5— function potential as a 
perturbation: 

n{^)=no{^)-vo6{^) (52) 
we consider the analogous differential problem: 



(53) 



G{L; X, y) = Go{L; x, y) " ^ ds J d'^zGo{L - s; x, 2)5{7)G{s; z, y) (55) 



G(0;x,y) = 5(x-y) (54) 

We wish to compute G{L] x, y) starting from the Green function Go{L; x, y), which is sup- 
posed to be known. It is possible to show that G{L] x, y) and Go{L; x, y) are related by the 
integral equation [sfll. l3l|: 

f-L 
10 

We see that in the right hand side of the above equation the presence of the 5— function 
forces us to consider the functions Go{L; x, y) and G{L; x, y) evaluated at the points x = 
and/or y = 0. Usually, at these points Green functions may be not well defined due to the 
presence of singularities. A concrete procedure to remove these singularities will be indicated 
later. For the moment, we go further with formal manipulations, assuming that some kind 
of consistent regularization of the possible divergences has been introduced. 
First of all, we perform the integration over (fz in Eq. (f55|l : 

Gm X, y) = GoiL; x, y) - Vo T dsGoiL - s; x, 0)Gis; 0, y) (56) 

JO 
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The integral in ds appearing in the right hand side of Eq. ()56|) is a convolution which can 
be better treated after a Laplace transform. Thus, we transform both sides of this equation 
with respect to L: 

G{E; X, y) = Go{E; x, y) - VoGo{E; x, 0)G{E; 0, y) (57) 

where 

r+oD 

G{E;^,y)= e-^^G(L;x,y)f/L (58) 



and 

Go{E- X, y) = / e-^^Go(L; x, y)dL (59) 
Jo 

At this point, it is easy to extract from Eq. (|F7|l the expression of G{E; x, y): 

r(T^ \ r(w \ Go{E;^,0)Go{E;0,y) 

G{E- X, y) = G,{E- x, y) - .l^^^^^.q^q) ^^^^ 

The above formula may be used in order to solve Eq. (|in|l . In this case, TIq coincides with 
the free action: 

H, = j-V' (61) 
and the function Go{E] x, y) is given by: 

Go{E- X, y) = -Ko{V2E~c\ic - y|) (62) 

TT 

Here Ko{z) denotes the modified Bessel function of the second kind of order zero. 

Clearly, we cannot apply directly Eq. ([6n|) without introducing a regularization. As a 
matter of fact, if not treated, the naive denominator in the second term of the right hand 
side is equal to infinity, i. e. ^ + Go{E]0,0) = +oo. This is due to the fact that Kq{z) 
diverges logarithmically in the limit z —>■ 0: 

Ko{z) log 2 for2~0 (63) 

A natural regularization is suggested by the fact that, in polymer physics, it has no sense 
to consider lengths which are smaller than the size of the molecules which compose the 
polymers. Thus, it seems reasonable to regulate ultraviolet divergences by introducing a 
cut-off a at short distances. The length a is comparable with the molecular size. According 
to this prescription, by inserting the Green function of Eq. (f62|l in Eq. (j6fl|l . we obtain: 
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The symbol = means that the quantity in the left hand side of an equation has been replaced 
in the right hand side with its regulated version. The above Green function is what we need 
in order to evaluate explicitly the amplitudes of Eqs. (j4TH43ll . 

The infinities coming from the 5— function potential should be treated with some care in 
order to avoid ambiguities. For this reason, we would like to compare the naive prescription 
used here to derive Eq. ([Bill with the more rigorous procedure of renormalization. It is 
known in fact that the renormalization of the infinities coming from 5— function interactions 
p.od..es p,..ca,„ se„.Ue .suits Q. The d.ve.,e„ces will be .eclated ..od.e.„, a 
cut-off A in the momentum space. As a consequence, it will be convenient to express the 
free Green function of Eq. (jH^ in momentum space. To this purpose, we use the following 
formula: 

X „ gjp-(x-y) 

Ko{m\^ -y\) = — £p (65) 



2tt J ' + 

To evaluate the Green function at the singular point x = y = we need to compute the 
following divergent integral: 



= - / (66) 



1 r d^p 
271 J + 

Using the above cut-off prescription to eliminate the ultraviolet singularities we get, in the 
assumption A^ m^: 

I(m) ~ log — (67) 

m 

Now, according to the spirit of renormalization, we subtract the infinities from the bare 
parameters of the theory. In our case, after choosing an arbitrary mass scale n, which gives 
the renormalization point, we renormalize the bare coupling constant Vq. Actually, it will 
be better to call it Vbare instead of vq in order to distinguish it from the effective coupling 
constant Vq appearing in Eq. ()64|l . The subtraction of infinities is performed in such a way 
that the quantity: 

Go i?;0,0 = + log — -—log — 68 

becomes finite. We choose a sort of minimal subtraction scheme, in which the renormalized 
coupling constant Vren is related to the bare coupling constant Vbare as follows: 
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Applying the last two above equations back to Eq. (f6n|l . we get as a result: 



Eqs. (f64l) and (f7fl|l are reciprocally compatible. In fact, since a is very small, because it 
is the smallest possible length scale in our polymer problem, one can use the following 
approximation (see Eq. (fM| ) in the denominator of the second term of Eq. 

- + -KoiV2E~ca) ~ - - 7^ log(2Eca) (71) 

Comparing with the analogous denominator in Eq. ((201), it is possible to relate a with the 
inverse of the mass /x: 

Moreover, the effective coupling constant vq of Eq. (j6i|l may be identified with the renor- 
malized coupling constant Vren, which gives the strength of the repulsive interaction Q at 
distance scales of order a. 

Before concluding this Section, we make a small digression about the translational invari- 
ance of the free Hamiltonian (f6T|) and consequently of the free Green function (|62|) . Clearly, 
this is not the same translational invariance that was already present in the original action 
(III) due to the translational invariance of the potential This new invariance is rather 
related to the fact that the physics of the two polymer system in the absence of any inter- 
action does not change if we modify the relative positions of the polymer ends at t = and 
t = L in a symmetric way. An example of such transformations is the translation of both 
ends of polymer A at the initial and final points by a constant vector a: 

r^(0) = TAiO) + a (73) 
TAiL) = rA(L)+a (74) 

As a result of the translations (jZSHZlI), the relative vector r{t) of Eq. ^ at the instants 
t = and t = L changes as follows: 

r'o = ro + a (75) 

r'l = n + a (76) 

Clearly the propagator ()62p is invariant under the above transformations. This kind of 
invariance can be explained as follows. As far as the two polymers A and B do not interact. 
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each of them may be treated as an independent system. If we translate for instance both 
ends of polymer A at t = and t = L in the symmetrical way shown by Eqs. (f73|) and (fzijl . 
the number of configurations of polymer A and consequently the configurational entropy of 
the whole system do not change, because the transformation is equivalent to a translation of 
polymer A in the space. Of course, this invariance disappears as soon as the two polymers 
start to interact or if they are entangled together. Indeed, if one adds to the free Hamiltonian 
(j6T|) a (5— function potential, the propagator stops to be translational invariant as shown by 
the Green function of Eq. (|6fl|) . which does not depend on the difference x — y. 



V. REPULSIVE FORCES AND WINDING ANGLES: QUALITATIVE AND 

QUANTITATIVE CONSIDERATIONS 

In principle we have at this point all the ingredients which are necessary to compute the 
second moment of Eq. (f^ . In Eqs. and fHnil . in fact, the quantities N{E;ri,rQ) and 
D{E; ri, tq) are written as linear combinations of the amplitudes of Eqs. (HTHiHll . which can 
be explicitly evaluated using the propagator G{E, u, v) given in Eq. j45j and the formu- 
las of Eqs. (|iil - li8|l . The remaining task is to perform the integrations over the coordinates 
X and y in Eqs. (|iil - li8|) . From the analytical point of view, the evaluation of these integrals 
poses severe technical problems, which can be solved only with the help of drastic approxi- 
mations. However, the difficulties become milder if we average the second moment over the 
endpoints of the polymers as shown in Eq. (fH^ . In the energy representation, which we 
are using, this means that we have to consider the following averaged version of the second 
moment in Eq. (f26j) : 



™ = urn p^' 



where 



N{E) = J J rfViiV(E;ri,ro) (78) 

D{E) = J d?r, J d'nD{E- r,,ro) (79) 

Accordingly, we need to integrate the quantities -ft'(ri,ro) and /t^(ri,ro), = 1, . . . , 4 of 
Eqs. (j33H3HI) with respect to ri and tq. Putting: 

K{E) = J t/Vo J rfViir(ri,ro) (80) 
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I^{E) = j dho j d\j^{vi,vo) cu = l,...,4 (81) 

we obtain from Eqs. ()i2|) and ()i3|) the following expressions of N[E) and D[E): 

4 

N{E) = 2K{E) - J2 UE) (82) 

UJ = 1 

D{E) = y'rfVorfViG(E;ri,ro) (83) 

It will also be convenient to split the propagator G{E; u, v) of Eq. into two contributions: 

G{E; u, v) = Go{E; u, v) + Gi{E; u, v) (84) 

where Go{E; u, v) is the free propagator of Eq. (jl^ . which is invariant with respect to the 
transformations (f75|l and (f76|l . while 



GiiE; u, v) = -\{E)KoiV2E^\u\)KoiV2E^\^^\) (85) 

In the above equation we have isolated in the expression of Gi{E; u, v) the factor: 

X(E) = --(- + -Ko{V2E~ca)) ' (86) 
vr Vt>o vr / 

It is clear that the origin of the term Gi{E; u, v) in the propagator is due to presence of the 
5— function interaction ^ in the polymer action (jl}. In fact, if fo = 0, this term vanishes 
identically. Thus, using the splitting of the propagator of Eq. |84ll . it is now possible to 
separate in Eqs. (jUHlHI) the contributions to entanglement given by repulsive forces. 

It seems natural to expand the quantities D{E), K{E) and Iu){E) defined in Eqs. (jZH]), 
(|8n|) and (|8T|l with respect to \{E) as follows: 

D[E) = D^''\E) + D^^\E) (87) 
K{E) = (E) + i^(^) (E) + ir(') (E) (88) 
UE) = li'\E) + I(J\E) + li'\E) + /(f)(E) (89) 

where the superscript (ra), with n = 0,1, 2, 3, denotes the order in X{E). There are no higher 
order terms with n > 4, so the above expansions are exact. 

It is easy to show how K{E) and the I^{Eys depend on the pseudo-energy E. After a 
rescaling of the integration variables ri,ro,x and y in Eqs. (f7H| and (fTHl) . one finds in fact 
that: 

ii-W (E) = A"(E)E-2ir(") n = 0, 1, 2 (90) 
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n = 0,1,2,3 



(91) 



where the factors K^^'^'s and the /^"-''s are functions of the parameters a and c, but not of E 
or Vq. In fact, the coupling constant Vq appears only inside the powers of X{E). Let us note 
in Eqs. fnOll and (pT|l the presence of the overall factor E~^ in Eqs. 1^01) and ^T^ . Looking 
at Eq. (|82|l . it is clear that the whole function N{E) is characterized by the leading scaling 
behavior N{E) ~ E~'^. In the L— space, after an inverse Laplace transform, this behavior 
corresponds to the following scaling law, which is typical of ideal polymers: N{L) ~ L. The 
powers of X{E), appearing in the expressions of K'^"'\E) and I^\E), introduce corrections 
to this leading behavior that are at most logarithmic rn. E. As a matter of fact, if the 
condition 2Eca? -C 1 is satisfied, we have that: 



Naively, the above seems the only logarithmic correction which is possible in the expressions 
of N{E) and D{E) when E is small. However, that this is not true. In fact, in deriving 
Eqs. (f90|) and (|9T|) . we have not considered the divergences which arise in some of the 
integrations over the variables x, y , Tq and ri . After regulating these divergences with some 
prescription, as for instance the ultraviolet cut-off a used in Eq. ((6Hl . we will see in Section lVll 
that the naive rescaling of variables exploited in order to obtain Eqs. (|nn|l and (^JTll does no 
longer work and one should add extra logarithmic factors to these equations. 

Eqs. (|90|l and (f9T|l may be also useful to study the case of polymers in confined geometries. 
As a matter of fact, for large values of E^ one recovers the limit of small values of L, in 
which the region between the initial and final height is very narrow. Looking at Eqs. ()90|) and 
(j9T| . it is clear that the only interesting corrections when E is large come from the powers 
of X{E). To evaluate these corrections, one should note that the modified Bessel function of 
the second kind Kq{z) goes very fast to zero for large values of z. As a consequence, already 
in the domain of parameters in which 2Eca? > 10, it is possible to make the very interesting 
approximation 



Unfortunately, it turns out that the values of the energy for which the above equation is 
satisfied are not physical, as it will be shown below. 

Other useful information on the influence of repulsive forces on the winding angle can 
be obtained studying the form of the function Gi(i?;u, v) of Eq. (|85|) . We remember in 




(92) 



X{E) ~ Vq 



(93) 



vr 
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fact that all the effects of the repulsive forces are concentrated in this component of the 
propagator. Supposing for example that the value of |u| is very large, i. e.: 



we have the following approximate expression of Gi{E] u, v): 



A relation analogous to Eq. may be written also for the variable v. In practice, Eq. 
means that the repulsive interactions do not play any particularly relevant role in polymer 
configurations in which the ends of the trajectories at some point are very distant. This 
is not a surprise. If the trajectories at some height t are very far from each other, they 
will have little or no chance to interact together via the repulsive interactions of Eq. (jSj), 
which are of short range. Eq. gives the concrete law with which the contributions of the 
repulsive forces are suppressed in configurations of this kind. In particular, if the distance 
between the trajectories is much greater than the characteristic length scale 



the infiuence of the repulsive forces ceases to be relevant. Of course, even if at some points 
the trajectories are very distant, polymers will always have a chance to get near enough 
to be able to interact if they are sufficiently long. As a consequence, we expect that the 
characteristic length Irep increases with the increasing of the lengths of the trajectories. It 
is easy to check that this is exactly the case. To show that, let us consider the dependence 
of Irep on the polymer length. One parameter which determines this length is the distance 
L between the ends of the polymers along the t— axis. Indeed, a trajectory connecting the 
two ends of a polymer must be very long if these ends are located at very distant heights. In 
the energy representation, large values of L correspond to small values of E. For example, 
in the limit E = 0, which corresponds to infinite polymer lengths, we have that Irep = oo, 
confirming our intuitive expectations. Another confirmation comes from Eq. (|1()4|) below, 
where a rough estimation of the behavior of Irep with respect to the distance L is given. The 
dependence on L is not the whole story. As a matter of fact, during their random walk in 
the t— direction from t = to t = L, polymers are also allowed to wander in the remaining 
two directions. Loosely speaking, the variations in the length of the trajectories associated 




(95) 




(96) 
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to the fluctuations in these transverse directions are taken into account by the parameter 

n 

c. Smaller values of c correspond to longer trajectories and vice- versa, see paf. It is now 
easy to realize from Eq. (f96| that the characteristic length Irep increases when c decreases 
as expected. Taking into account all the above considerations, it is possible to conclude 
that repulsive forces give relevant contributions to the second moment only in the case of 
configurations of the system in which the trajectories of the two polymers are not too far 
from each other. As a matter of fact, in the propagators appearing in the amplitudes of 
Eqs. (I33HiHI) all configurations in which the distance between the trajectories at some height 
in the t— axis is bigger than a few characteristic lengths Irep are exponentially suppressed 
according to Eq. One may also add that this suppression becomes milder in the case 

of long polymers, because we have seen that the characteristic lengths Irep grows with the 
length of the polymers with a law which has been given in Eq. ()in4p . 

In the rest of this Section we will analyze some interesting limiting cases, in which repul- 
sive interactions become particularly weak or strong. To this purpose, it would be appealing 
to consider the quantity c~^X{E), where \{E) has been given in Eq. as an energy 

dependent effective or running coupling constant of the repulsive interactions. This could 
be suggested by the expansions of Eqs. (|871 - l89| and by the fact that the quantity c^^X{E) 
has the right dimension to be a coupling constant. Indeed, we will see that there are cases in 
which the strength of X{E) really determines the strength of the repulsive forces. However, 
this is not true in general, as it should be because \{E) is just a parameter which has been 
factored out from the expression of Gi{E; u, v) and thus its meaning does not coincide with 
that of a running coupling constant. Keeping that in mind, we start to study the perturba- 
tive regime, in which Vo is very small. In the part of the propagator in which there are the 
contributions of the repulsive forces, i. e. the function G'i(i?;u, v), vq is present only inside 
X{E). Expanding this quantity in powers of vq, we obtain: 



We see that, at the leading order in Vq, X{E) is proportional to Vq and thus, as it could 
have been expected, Gi{E; u, v) may be treated as a small perturbation with respect to the 
free propagator Go{E,u,\). Let's now go back to Eq. In that equation it turns out 

that X{E) has the same behavior as in the perturbative regime, even if Eq has been 
derived in the hypothesis that 2Eca? > 10, but without supposing that vq is small. Before 




(97) 
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dwelling upon the physical meaning of this coincidence, let's see what is the significance of 
the condition 2Eca? > 10. To this purpose, we make the following approximations: 

E-^ - ~ a (98) 

c 

As mentioned before, it is quite reasonable to assume that the length L is proportional to 
the inverse of the energy E, while the second approximation implies that polymers are very 
flexible. For example, in polyethylene the Kuhn length~ 1/c is of the order of molecular 
sizes. Exploiting Eq. (f98|l . it turns out that the condition 2Eca? > 10 is equivalent to the 
condition L < ^. This would mean that our system is squeezed in a volume whose height L 
is smaller than the size of a monomer. Clearly, this situation is not very physical. 

Since we have been able to compute the exact form of the propagator G{E; u, v), it is not 
difficult to study also the strong coupling limit Vq — > oo. As in the perturbative case, the 
only affected part of the propagator (j6i|l is the factor X{E) appearing in Gi{E; u, v). After 
a trivial calculation, one finds that, in the strong coupling limit, the form of Gi{E] u, v) is 
given by: 

G'i(E;u,v) ~ (^Ko{V2E^a)y^ KoiV2E^\u\)Ko{V2E^\^r\) (99) 

Assuming that polymers are very long, let us study the left hand side of the above equation. 
This is a ratio of modified Bessel functions of the second kind. Since a is a very small 
quantity and these functions have a logarithmic singularity if their argument is small, see 
Eq. it is licit to suppose that 

Ko{V2E^a) > Ko{V2E^\u\)Ko{V2E^\^r\) (100) 

unless |u| ~ a and/or |v| ~ a. On the other side, we know from Eq. (fQijl that, if 
|u|, |v| ^ the product of modified Bessel functions Kq{^/2Ec\u\)Kq{\/2Ec\y\) de- 

cays exponentially. In other words, in the left hand side of Eq. (j99l) the denominator will 
dominate over the numerator whenever the distance between the polymer trajectories is not 
of the order of a few molecular sizes. Thus, if fo is large, the major contributions to winding 
angle coming from the repulsive interactions occur when the trajectories are very near to 
each other. This could be expected from the fact that, in the strong coupling limit, one 
recovers the excluded volume interactions. 

Finally, let us study the domain of the parameters E and c in which the condition 

2Eca^ < 1 (101) 
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is verified. We will see that this domain is particularly interesting, because if condition ()lflip 
is verified, the corrections of the repulsive interactions to the entropy dominated behavior 
of ideal polymers become relevant. It has been already shown that under the assumption 
made in Eq. (IIOIL the parameter X{E) is approximated as in Eq. (jH^ . Even if it is not 
strictly necessary, we suppose here that vq has some finite value, while polymers are so long 
that the following inequality is satisfied: 

< -- log (V2E^a) (102) 

This further assumption is to eliminate the dependence on fo, which could introduce con- 
fusion in the following discussion due to possible interferences of condition ()inij) with those 
of the perturbative and strong coupling regimes. In the L— space, Eq. (jin2jl corresponds to 
the inequality e^'^l''^^ -c ^. Now G\{E\ u, v) may be approximated as follows: 

G^{E- u, v) ~ / iro(v^|u|)iro(\/2^|v|) (103) 
vr log(V2-E ca) 

As promised, the above equation does not contain the parameter vq. We see from the left 
hand side of Eq. (|1()3|) that the function Gi(_E;u, v) is logarithmically suppressed, due to 
the presence of log(\/ZEca) in the denominator. This suppression effect is counterbalanced 
only at short distances by the two modified Bessel functions of the second kind appearing in 
the numerator, which diverge logarithmically whenever V2-E'c|u| = and/or \plEc\w\ = 0. 
The total result of these opposite effects in the expression of the averaged second moment 
will be presented in Section ED after performing the explicit computation of the amplitudes 
of Eqs. JSHIHI). 

To conclude this Section, let us give some concrete values of the involved parameters. 
First of all, let us estimate the values of L, for which the two polymer system is in the 
regime pOlll . Using the approximations made in Eq. (j98|) . we may conclude that, if the 
relation ()1()1|] is satisfied, the length L needs such that L ^ 2a, i. e. L is at least of the 
order of hundred molecular lengths or more: L > 100a. Moreover, it is possible to give a 
rough estimation of the maximum distance of the end points, after which the two polymers 
are too far from each other to allow a relevant contribution to the winding angle due to 
repulsive interactions. Using Eq. (|9i|) . in fact, it turns out that the repulsive interactions 
are relevant only in the range of distances: 




(104) 
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Finally, the situation opposite to condition ()lfllj) is not realistic, because it leads to the 
constraint L <^ 2a. This would corresponds to the case of a polymer which is shorter than 
the size of the molecules composing it. 

VI. CALCULATION OF THE AVERAGED SECOND MOMENT 

At this point we are ready to compute the quantities N{E) and D{E) of Eqs. (TZSl and 
(17^ . We start with D{E). Using Eqs. ^ and the splitting ^ of the propagator, 

one has at the zeroth order in X{E): 

D^^\E) = I d^o I rfViGo(i?;ri,ro) = J d^o j d^n^K^{V2E'c\v, - (105) 

After a shift of variables, the above equation gives: 

= S ! d^n-Ko{V2E^\ri\) (106) 

where S = J d^r^ is the total surface of the system in the two dimensional space, which is 
transverse to the t— axis. Using the identity 

/ rfVi-iro(v^lril) = 4 (107) 

one finds: 

D(°)(E) = S/E (108) 

This expression of D^'^\E) has the following interpretation: We are performing here an 
average of the second moment with respect to all possible initial and final positions of 
the endpoints of the polymers and D{E) counts the number of these configurations. The 
component D^^\E) of D(E) depends only on the free propagator Go{E;x,y), which is 
translational invariant in the sense discussed after Eq. (j62|) . This invariance explains why 
the number of configurations grows proportionally to the surface S. The reason is that, for 
each configuration of the polymers, one can obtain other equivalently probable configurations 
by the symmetric translation of their ends on the surface S at the initial and final instants. 
Let us now apply to D^°\E) an inverse Laplace transform, in order to go back to the 
L— space. After a simple calculation we obtain: 

D^''\L)=S (109) 
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i. e. D^^\L) does not depend on L. 

The next and last contribution to D{E) is given by: 

D('\E) = I dVorfVG'i(i?;ri,ro) = / £rod'n^X{E)Ko{V2E~c\r,\)Ko{V2E~c\ro\) (110) 
Exploiting Eq. (|in7|l to integrate out the variables Tq and ri, we get: 



D^'\E) = -\{E)E 



(111) 



We remark that the above contribution to D{E) vanishes in the limit fo = 0. This could be 
expected due to the fact that D^^^E) collects all contributions coming from the repulsive 
interactions. These interactions break explicitly the translational invariance of the free 
part of the action and, as a consequence, D^^\E) is no longer proportional to the surface 
S as D^^\E). Unfortunately, it is not easy to compute the inverse Laplace transform of 
D^^\E) without making some approximation. To this purpose, we assume that the repulsive 
interactions are weak, i. e. fo <^ 1, and that the value of L is large. In this case, since we 
are in the domain of small i?'s, it is possible to expand D^^\E) up to the second order in 
Vq as follows: 

vr /c_. fc , ^ ^ ^^^^^ 



nm^E)r.-{-E 
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C \ TT 



^^0 



C 

-Vo 
TT 



E-^log{V2E 



ca] 



In order to obtain the above equation we have used both Eqs. (jH!^ and (jHIj). The inverse 
Laplace transform of Eq. ()112|) gives: 

C-V 



TT 



^0 "^0 log(v2ca) + 



L + — u^LlogL 



(113) 



where C ~ 0.577215664 is the Euler constant. 

Putting Eqs. (jl()8|] and (llllll together, we obtain: 



D{E) = D^''\E) + D^^\E) = SE-^ + -X{E)E 



(114) 



This is an exact result. An approximated expression of D{L) can be derived instead from 
Eqs. (ITOQl) and dm . 

Now we turn to the derivation of N{E). We start by computing order by order in A(-E') 
the contributions to the quantities K{E) and Iuj{E) of Eqs. (|HHll and (|Hn|l respectively. At 
the zeroth order we have for K{E): 

K^o\E) = ^J rf'xA'(x) J (friKo{V2E'c\Ti - x|) ^ d\oKo{V2E^\ro - x|) (115) 
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After performing an easy integrations over the coordinates ro,ri, one obtains: 



2c 

The remaining integral with respect to the x coordinate is both ultraviolet and infrared 
divergent and needs to be regulated. We have already seen that the singularities in the 
ultraviolet domain may consistently be eliminated with the introduction of the small dis- 
tance cut-off a. A large distance cut-off is instead motivated by the fact that the size of 
a real system is necessarily finite. Implicitly, we have already used this kind of infrared 
regularization in Eq. where we have assumed that the total surface S of the system 

in the directions which are transverse to the t— axis is finite. Supposing that the shape of S 
is approximately a disk of radius R, so that S ~ irR'^, we may write: 



/ 



d^xA\^) = - -f (117) 



Substituting Eq. UllH in Eq. (I116|) . one obtains the following expression of K^^\E): 



K^'\E) = ^E'Hog(^) (118) 
8nc Ka'^n ' 



The inverse Laplace transform of K'^'^^E) gives: 



ir(°)(L) = AlogfA) (119) 



6TTC ' Va^TT 



We have now to compute the quantities /^°)(£^), with cu = 1, ... ,4. The expressions of 
the /^°^(-E)'s may be obtained from Eqs. (P?T| and pKHi^ . by substituting everywhere the 
propagator G{E; x, y) with its free version Go{E; x, y). It is easy to show that: 

li°\E) = forcj = 1,...,4 (120) 

This vanishing, which is actually a double vanishing, is due to the fact that each of the 
/^°''(-E)'s contains an integral of a total divergence together with an integral which is zero 
for symmetry reasons. For some values of u, like for instance when u = 3, to isolate such 
integrals it is necessary to perform some integrations by parts. This is allowed because the 
/^°^(-E')'s are not affected by divergences, contrarily to K(E). 

As an example, we work out explicitly the case of Ii'\e). The first vanishing integral is 
the following: 

J d'roViGoiE; ro, y) = ^ / rfVoVjiro(\/2^|ro - y|) (121) 



28 



This is of course zero due to symmetry reasons. The second vanishing integral in Ii'\e) is 
of the form: 

1 = Jd^xJ d\iA{^)Go{E; n, x) V5,Go(^; y, x) (122) 

After performing the integration over ri with the help of a shift of variables and of Eq. (110711 . 
we have, apart from a proportionality factor: 

/ oc y"d2xAi(x)V^Go(^;y,x) (123) 

Since v4j(x) is a divergenceless vector potential, i. e. Vx^i(x) = 0, / can be rewritten as the 
integral of a total divergence: 

1=^1 d\\/l, {Ai^)KoiV2E~c\y - x|)) (124) 

Clearly, the left hand side of the above equation is zero. This fact can be also checked 
passing to the Fourier representation. Exploiting Eq. (f65|) and the formula 

in Eq. (jl23jl . one obtains: 

Thus / = because e^^piPj = 0. In an analogous way one shows that also l2'\l^^^ and I^^ 
are identically equal to zero. 

We are now ready to compute the contributions to N{E), which are linear in \{E). First 
of all, we treat the term K^^\E), which is given by: 

K^'\E) = ^J d^xj d^ro I dViA2(x) [Gi(i?;ri,x)G'o(i?;x,ro) + G'o(i?;ri,x)Gi(i?;x,ro)] 

(127) 

The integrations over tq and ri may be easily performed using Eq. (Il()7|] and give as a result 
a factor which is proportional to E'"^. After that, only the following integral in x remains 
to be done: 

/ dW(x)Ko(v^|x|) = / d^x-^KoiV2E~c\^\) (128) 

J l^^J •^|x|>a |x| 

Here the ultraviolet divergence, which is present in the left hand side, has been regulated 
in the usual way with the introduction of the short distances cut-off a. Infrared divergences 
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are absent. Going to polar coordinates, the right hand side of the above equation becomes: 

' d^x—K^{\f2Ed(s\) = — / dp^^ ^ (129) 



(27r)2 7|x|>a |X|2 ' 2lT. 'a p 

Putting everything together, one arrives at the final result: 



= J-E-H(E) r (130) 

2,7CC Ja P 



If the quantity y2Eca is small, it is possible to derive the following asymptotic expression 
of ir«(E): 

(E) ~ —E-^\{E) log2( V2^a) (131) 
Attc 

To go from Eq. ()13np to Eq. (|13H1 . we have used the asymptotic formula: 



+00 



Ko{V2Edp) 1 



a 



dp^^ ^ ~ -log2(V2E^a) (132) 

P ^ 



which is valid for small values of y/2Eca. We see from Eqs. (I13()jl and (jl32|] that the 
presence of ultraviolet divergences, together with the needed regularization, has modified 
the naive form of K^^^E) function of the pseudo-energy E given in Eq. ()9fl|l . The 
modification consists in the appearance of the factor /^°° —Kq{\/2Ecp), which exhibits a 



square logarithmic singularity in the limit \f2Eca = 0. 

The inverse Laplace transformed of K^^^ (E) can be derived only making some approxi- 
mation. As in the case of D^^\E)^ we will work in the double limit, in which vq is very small 
and L is very large. After a few calculations we obtain: 

K^'^ (L) - £ ds [log(L -s) + C] (log s + C) + 

jL\og\2ca^) + ^log(2ca2) [(C - 1)L - LlogL]| (133) 

At this point we have to compute the expressions of the /^^^(i?)'s, cu = 1, . . . ,A. It is 
possible to show that these contributions vanish identically, i. e.: 

jW(E) = forcj = 1,...,4 (134) 

The motivations of this vanishing are similar to the motivations for which there are no 
contributions at the zeroth order: All terms which appear in the quantities I^\E) contain 
at least one integral of a total divergence or one integral, which is zero for dimensional 
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reasons. As in the case of the /j^°)(-E')'s, there are some values of u for which it is necessary 
to perform an integration by parts in order to isolate these vanishing integrals. Once again, 
this is allowed because the II^\E)^s do not contain divergences. 
At the next order in A(£'), we have the last contribution to K{E): 

K^^\E) = 1-Jd'xj d'nj dVoA2(x)G'i(E;ri,x)Gi(E;x,ro) (135) 

After performing the integrations in ri and Vq with the help of Eq. ()107p . Eq. ()135p becomes: 

K^^){E) = ^X\E)E'^ J rfW(x) (Ko(v^|x|))' (136) 

The integral in x is divergent and needs a regularization. Going to polar coordinates, we 
obtain the result: 

K(^)^E) ^ -\'{E)E-' ^ (Ko{V2E~cp)Y (137) 

C J a p ^ ' 

Also in this case, we note that the presence of the regularization modifies the dependence 
of K^'^\E) on the pseudo-energy E with respect to the naive formula of Eq. (jHOI)- The 
correction consists in the factor J^°° ^ (^Ko{y/2Ecp)^ . In the limit \/2Eca = 0, this factor 
diverges as powers of \og{\/2Eca). 

To conclude the analysis of the contribution to N{E) at the second order in \{E), we 
show that the II^\E)^s are identically equal to zero. As a matter of fact, it is easy to verify 
that for ui = 1, 2, 4 each I'^\E) contains terms of the kind: 



5(x) = A(x)V^i^o(v2Ec5|x|) (138) 
These terms vanish identically because of the following identity: 



1 


X' 


' dKa{^/2Ec\ 


|x|) 


^/2Ec\ 


x 


|2 9|x| 





V^^K,{V2Ec\^\) = . r (139) 



Substituting Eq. (11391) in Eq. (I138|] and using the explicit expression of the vector potential 
A(x) of Eq. ((H), we get: 



1 ejix'x^ dKo{^/2Ec\ 


|x|) 


2'kV2Ec 


x 


|4 9|x| 





Clearly, the right hand side of the above equation is zero because ejiX^x^ = 0. If c<j = 3, 
instead, the vanishing function -B(x) of Eq. ()138|) may be isolated in the expression of 
lf^{E) = only after an integration by parts. 
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Finally, at the third order in X{E) we have only the quantities /^^^(-E')'s, since K{E) has 
at most quadratic powers of X{E). It is easy to realize that: 

j:iL'\E) = (141) 

UJ=1 

because the following relations hold 

I?\E) = -I^iE) = Jf (E) = -li'\E) (142) 

As a consequence of Eq. (jHlI), it is clear that there are no contributions to N{E) at this 
order. 

Using Eqs. (TTra . (THTO and jnH), we arrive at the final result for N{E): 

N{E) = -LE-Hog(^) + —\iE)E-^ ^KoiV2E~cp) 
Attc Va vr/ ttc J a p 

+ ^X'iE)E~' ^ (KoiV2E~cp)Y (143) 

C Ja p ^ ' 

We can now insert in the formula of the second moment of Eq. (|77|l the functions -D(i?) 
and iV(-E) given in Eqs. (jll4jl and (jl43jl respectively. The outcome is: 



E- 

{m'){E) = - 



4^ log (J,) + A(E)i:^::i^^ + x^e^-^ jr ^ {mvm-cp) 



S + l\{E)E-^ 

(144) 

In the L— space, the already mentioned difficulties with the computation of the inverse 
Laplace transform of D{E) and N{E) allow an analytical result only in the double limit of 
weak coupling constant Vq and of large values of L. At the first order in fo, the expression 
of {rn?) reads as follows: 

^ = S^L ^'^'^ 



where K^^\L) has been given in Eq. jn 

So far, we have considered the averaged second moment of Eq. (f26|l . corresponding to the 
case in which the polymer ends are not fixed. In the energy representation, we have seen 
that this version of the second moment can be exactly computed. To conclude this Section, 
we would like to show that it is possible to provide also an exact expression of the second 
moment {'m'^)ri,ro in the L— space and with fixed polymer ends up to an inverse Laplace 
transform of the propagator given in Eq. ((Hlj). The starting point is the exact formula of 
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the second moment (m^)ri,ro(-E') of Eq. (|i9|l . All the ingredients of this formula are defined 
in Eqs. and Looking at Eq. (gni), it is clear that: 

4 

N{E- n, ro) = 2K{ti, tq) - ^ /.(ri, ro) (146) 

and 

D{E;ri,ro) = GiE;ri,ro) (147) 

Let us note that the functions /a;(ri, ro) are all equal up to integrations by parts, which can 
shift the differential operators Vx and Vy in Eqs. (j44ll - (|48ll . This fact will be used in order 
to simplify the expression of the inverse Laplace transformed of N{E] ri, tq) in the L— space. 
To compute the inverse Laplace transforms of both A^(-E'; ri, ro) and D{E; ri, ro), we use the 
following property of the inverse Laplace transform of the product of two functions f{E) 
and g{E): 

C-\f{E)g{E)) = f^dsfiL - s)g{s) (148) 
Jo 



Applying Eq. ()148|] to evaluate the inverse Laplace transforms of K{E) and of the /(^(ri, ro) 
in Eqs. (I146|) and (I147|) . we obtain after some calculations: 

N{L;ri,ro) = - (i^xA^(x) / (isG'(L — s; ri, x)G(s; x, ro) 
c J Jo 

- ^ J (fx J (fy J G'(L-s;x, ri) 

X l^ds'didlG{s-s';y,^)G{s';ro,y)A{^)My) (149) 
D(L;ri,ro) = G(L;ri,ro) (150) 

The second term in the right hand side of Eq. ()149p is the contribution given by the func- 
tions /t^(ri,ro), uj = 1, ... ,4, while the first term comes from -ft'(ri,ro). Remembering the 
definition (f2T|l of the second moment in terms of iV(L; ri, ro) and D{L] ri, ro), we get: 



(m2)r„ro = [G(L; n, ro)]"^ 



9 r fL 

C 



J (i^xA^(x) J dsG{L — s;ri,x.)G{s;:x.,ro) 



^^Jd'x J d'y dsG{L-s;^,r,)J^ds'dldlG{s-s';y,^)G{s';vo,y)A{^)My) 



(151) 



If we knew how to compute the propagator G{L; x, y) starting from its Laplace transformed 
(jnH), we could evaluate explicitly the expression of the second moment in the L— space. 
Unfortunately, it is too complicated to perform the inverse Laplace transform of the propa- 
gator G{E]x,y). Due to this technical difficulty, Eq. ()15ip is only formal. Progress can be 



33 



made however in the limit fo = 0, in which the propagator is given by the Green function 
Go{E] X, y) of Eq. (f62|) . This will be done in the next Section. 

VII. THE CASE OF IDEAL POLYMERS 

In order to allow the comparison with previous results, this Section is dedicated to the 
case of ideal chains in which Vq = 0. First of all, we discuss the formula of the averaged 
second moment derived in the previous Section, Eq. ()145p . In the limit fo = 0, Eq. ()145p 
becomes: 

The presence of a geometrical factor like the surface S of the system in the expression 
of (m^)o has been already related to the translational symmetry of Eqs. ((7511 and (f76ll . 
Assuming that this surface has approximately the shape of a disc of radius R, we can put 
S = ttR'^ as in Eq. (IllTj) . Eq. (jl52|) predicts that the average degree of entanglement scales 
as follows with respect to the distance R: 

[m )o oc (153) 

The meaning of Eq. ()153|1 is the following. We remember that the averaged second moment 
(m^)o describes the entanglement of two closed polymers whose ends on the surfaces at 
t = and t = L are not fixed. In this way, the polymers are allowed to move freely and it is 
natural to suppose that, the bigger will be the volume SL in which the polymers fluctuate, 
the bigger will be also the average distance between them. Thus, if the surface S increases 
its area, the probability of entanglement must decrease. The exact law of this decreasing is 
given by Eq. ()153p . 

On the other side, one would expect that the probability of getting entangled is higher 
for long polymers than for short polymers. Eq. ()152ll gives a result which is in agreement 
with the above expectation, because the second moment (m^)o scales as follows with respect 
to the parameters L and c, which determine the polymer length: 

{m^)o oc - (154) 
c 

In particular, one can show that the total length of a polymer increases proportionally to 
L and it is inversely proportional to the square root of c j^^. Accordingly, we see from 
Eq. ()154p that (m^)o increases proportionally to L and inversely proportional to c. 
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At this point we wish to study the second moment (m^)o,ri,ro of polymers with fixed 
endpoints. The subscript has been added to the symbol of the second moment to remember 
that we are working in the limit Vo = 0. Since we are dealing with ideal polymers, we have 
to substitute everywhere in Eq. the full propagator G{L; x, y) with the free one. The 

result of this operation is: 



2 f 

- / (fxA^hi) / dsGol-^ — ■s; ri, x)Go(s; X, ro)- 
c J Jo 



^Jd^xJ d% dsGo{L - s;^,ri) ds'didiGois - s';y,^)Go{s']ro,y)A{^)A,{y) 



(155) 

We notice that, as it could be expected, Eq. 1)155^ coincides with the expression obtained in 
for the second moment of one polymer winding up around an infinitely long straight wire 
lying along the axis. Luckily, the propagator Go{L;ri,ro) can be explicitly constructed 
upon computing the inverse Laplace transform of the propagator Go{E; ri; tq) of Eq. (f62ll : 

G'o(I^;x,y) = ^e^l--^l' (156) 



It is easy to check that the second term in the right hand side of Eq. ()155jl . which is 
associated with the contributions coming from the Iuj{ri, ro)'s, does not grow with increasing 
values of L. As a matter of fact, after a rescaling of variables, the numerator of this term 
gives: 

2 r r rL 



C2 



J (fx J d^y dsGoiL - s;x,ri) dsdld^Gois - s';y,x)Go{s';ro,y)Ai{x)Aj{y) 



An^L J J ^ Jo 1-t Jo t-f 



J- e U,/ '•Q P 

-e 2(^1^ ^ I (157) 

In the limit L — > oo, the quantity in the right hand side of the above equation scales as 
AL^-^, where A is a constant. Moreover, the propagator ()156|1 . which is in the denominator, 
scales as L^^. Thus, the ratio between the right hand side of Eq. (jl57jl and the propagator 
(115611 does not depend on L. This completes the proof of our statement. 

As a consequence of this statement, as far as the scaling of (m^)o,ri,ro for large values of 
L is concerned, it is possible to make the following approximation: 

("^^)o,ri,ro ~ -[G'o(-^^;ri,ro)]"W (i^a;A^(x) / dsGo{L - s;ri,x)Go{s;:ic,ro) (158) 
c J Jo 
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Unfortunately, despite of the fact that we are treating ideal polymers, the integral in (Px 
appearing in the above equation is still complicated and requires some approximation to be 
evaluated analytically. We will apply to this purpose the strategy used in Ref. Q to compute 
the second moment of three dimensional polymers, adapting it to our two-dimensional case. 
First of all, let us note that the integral in (I158|] is ultraviolet divergent. However, the 
infrared divergences which appeared in the energy representation are absent. This is due 
to the behavior of the propagator Go{L;:x.,y), which is much milder at infinity than the 
behavior of the Green function Go{E; x, y). To regulate the singularities at small distances, 
we proceed as usual by introducing the cut-off a. After a rescaling of all variables similar to 
that of Eq. (ITFtII . we get: 

2 

("^^)o,rl,ro ~ -[Gol^; ri, ro)]"^ x 

(Px' 1 ds' / C \2 1 |x'~ri,/^|^ _i|x'-rn /^f 



\x'\>^ L x'2 Jo s'{l - s') \2tt 
To go further, following Q], we assume that the relevant contribution to the integral in d'^x' 
comes from a narrow region around the singularity in x' = 0. Thus, we may put 



d'^x' L 



g 2{l-s''7l'^ '"^Vil g 2{l-s') 




27rlog \\l-a\e W^)\''^V^\ e W^l^^V^l (160) 



After making the above crude approximation, we obtain: 

g 

("^^)o,ri,ro r[G'o(^; ri,ro)]~Mog 

TlL 

1 




— e 2(1-^^') ^L^-J^^OL 



.g 2(l~s') ILg 2s' ^)L-)- 

(161) 



In deriving the above equation we have used the simple relation = (Yz^tj + 7- Let us 

now study the integral 

1 ds' 



1= / — e 2(i-.')'iie-2F'-0L (162) 
Jo s' 

The other integral in ds' appearing in ()16HI can be treated in the same way after the change 
of variables 1 — s' = t. It is not to allowed to take in the right hand side of Eq. (I162|] the 
limit L — > 00 because in this way the integral will not be convergent due to the singularity 
in s' = 0. For this reason, we split the domain of integration as follows: 
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+ / _e W^)^^L^-^^oT (163) 

Ju S 

where < m < 1. Clearly, the second integral converges after performing the limit L — > oo 
in the integrand and gives: 

The first integral instead diverges logarithmically with growing values of L. However, now 



it is possible to expand the exponential e il Jj^ powers of its argument, because the 

singularity in s' = 1 lies outside the interval [0, u]. Keeping only the leading order term with 
respect to L, we get: 

~ / UY^ c\ 

/~-Eif-^j -logM (165) 

where Ei(2;) is the exponential-integral function. When L is large, this function may be 
approximated as follows: Ei(2;) ~ log{—z) and, as a consequence: 

/~-log(^) (166) 

The second integral which we have left in Eq. ()16H1 gives the same result. Putting everything 
together in the expression of the second moment of Eq. (jl61jl . we obtain the final result: 

(m2)o,r„ro ~ -2 log [^<^ log ~ 2 (logL)' (167) 

This is exactly the behavior of the second moment derived in Ref. 



VIII. CONCLUSIONS 



In this article we have studied the entanglement of two directed polymers from a non- 
perturbative point of view. Our formulas of the second moment, a quantity which measures 
the degree of entanglement of the two polymers, take into account the repulsive forces acting 
on the segments of the polymers and are exact. The averaged second moment defined in 
Eq. IIHI), a version of the second moment corresponding to the situation in which the end 
points of the polymers are free to move, has been computed in Eq. (I144|) as a function of the 
chemical potential E conjugated to the distance L between the end points in the t— direction. 
The case of free ends is relevant in the treatment of nematic polymers and polymers in a 
nematic solvent 0|. Let us note that also the expression of the second moment without any 
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averaging and in the L space can be computed. This has been done in Eq. ()15ip . However, 
this equation is exphcit only up to the inverse Laplace transform of the propagator (jMI), 
which is too hard to be obtained in closed form. 

Eq. (I144jl shows that the averaged second moment is of the form {m?){E) = E^^f{E). 
The overall factor coincides with the scaling power law of two ideal polymers. The cor- 
rection f{E) to this fundamental behavior due to the repulsive interactions is a complicated 
function of E, whose analysis would require numerical methods. Nevertheless, it is possible 
to identify a dominance of the repulsive interactions in the domain of parameters in which 
the condition \j2Eca ~ is satisfied. This corresponds roughly speaking to the situation in 
which polymers are very long. In this region, the scaling laws with respect to the energy E of 
the numerator and denominator appearing in the right hand side of Eq. (jl44jl are corrected 
by factors which are logarithmic powers of log(v2£^a), see for instance Eq. ()132p . 

One advantage of our approach is that it is easy to separate within the expression of the 
second moment the contributions of purely entropic origin which are typical of free polymers 
from the contributions coming from the presence of the 5— function potential in the polymer 
action. This is essentially due to the splitting ([HKI of the propagator G{E\ u, v) appearing 
in the amplitudes (|iil - li8|) . The component Gq{E\ u, v) of the propagator coincides with the 
propagator of ideal polymers, while the component G\{E\ u, v) takes into account the effects 
of the interactions. Thanks to the splitting (|95|) . it has been possible to study the way in 
which the repulsive forces affect the average degree of entanglement of the two polymers. 
This has been done in Section El Our results are in agreement with the intuition. The 
precise law with which the effects of the repulsive forces on the entanglement decrease when 
the distance between the trajectories increases is given by Eq. (j95|l . In Section El it has 
been discussed also the strong coupling limit, which should be taken to recover the limit 
of excluded volume interactions. In our exact approach, it is not difficult to consider the 
case in which the coupling constant fo is large. For instance, the component Gi(-E';u, v) 
of the propagator, which is responsible of the effects due to the repulsive interactions, has 
been given in the strong coupling limit in Eq. (^^1) . Studying the form of this component 
assuming that polymers are very long, it has been argued that, at strong coupling, the major 
contributions to the winding angle coming from the repulsive interactions occur when the 
trajectories are very near to each other. Many other qualitative and quantitative character- 
istics of the behavior of the two polymer system under consideration have been presented in 
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Section El 

The case of ideal polymers, in which fo = 0, has been discussed at the end of Section lVIIl in 
order to make comparison with previous works. The scaling of the averaged second_moment 
for large values of L obtained in Eq. (I154jl is in agreement with the results of |2], if one 
takes into account the fact that, after the averaging procedure of Eq. (fll?! and the infrared 
regularization of Eqs. (|in6|l and ()117|1 . one is effectively treating a system of polymers 
confined in a cylinder of finite volume SL. In Section IVIII we have evaluated the second 
moment, always of two ideal polymers, using the approach of Ref. j^. The outcome of this 
calculation, namely the scaling behavior of (m^)o,ri,ro ^it the leading order in L, is reported 
in Eq. (11671) . This result is in agreement with the square logarithmic behavior obtained in 

Q Q 

but not with the logarithmic behavior predicted in [7|. However, this discrepancy can 
be expected due to the fact that, in Section IVTH we have assumed, following Ref. that 
the most relevant contribution to the second moment coming from the integral in Eq. ()159p 
is concentrated in a narrow region near the singularity in x' = 0. This clashes with the 

n 

assumptions of Ref. |2|, in which instead it is argued that the main increase in the winding 
angle does not occur when the polymer trajectories are near, but rather when they are far 
one from the other. Finally, there is also an apparent discrepancy between the linear scaling 
with respect to L of the averaged second moment {m?)^ and the square logarithmic scaling 
of the second moment (m^)o,ri,ro- This disagreement is explained by the fact that, in the 
first case, the ends of the polymers are free to fiuctuate, while in the second case they are 
fixed. It is therefore licit to expect that two polymers with free ends are more likely to 
entangle than two polymers whose ends are constrained. 

Concluding, we would like to discuss possible further developments of this work, together 
with some problems which are still left open. First of all, the number of entangling polymers 
has been limited to two. To go beyond this restriction, one should explore the possibility of 
replacing the external vector potential Aj(x) of Eq. ^ with Chern-Simons fields. Abelian 
Chern-Simons field theories have been already successfully applied in order to impose topo- 
logical constraints to the trajectories of an arbitrary number of closed polymer rings in |37l |. 
We hope to extend those results also to the case of directed polymers in a forthcoming pub- 
lication. Of course, if the polymer trajectories are open, the constraints among them are no 
longer of topological nature as in [^^li so that the application of Chern-Simons field theory 
to directed polymers should be considered with some care. 
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We have also not made any attempt to introduce in the treatment of polymer entangle- 
ment more sophisticated constraints than those which can be imposed with the help of the 
winding angle. This is in effect still an unsolved problem, desp ite the fact that two pow- 
erful and strategies have been proposed for its solution ya, ISSlkOl- In the first approach, 
pioneered independently by Kleinert, Kholodenko and one of the authors [121, 123, 122l| • the 
constraints are expressed via the Wilson loop amplitudes of non-abelian Chern-Simons field 
theories. Some progresses toward a concrete realization of this program in polymer physics 
have been made in Refs. Illl,!^!- In the second approach, developed by Nechaev and cowork- 
ers, see |4fl,] and references therein, polymer trajectories are mapped on a complex plane 
with punctures. The link invariants necessary to impose the constraints are then constructed 
using the properties of conformal maps. 

Another possible development is the treatment of attractive interactions, in which the 
strength Vq in Eq. Q takes negative values. In this case, the 5— function potential admits a 
bound state 0| and the propagator of Eq. (j64ll develops a singularity, in which \{E) = oo, 
at the energy corresponding to this bound state. It would be extremely interesting to 
investigate how these facts affect the polymers' entanglement. Another issue which deserves 
attention is that of hairpin turns. Hairpins are important in nematic solvents and can 
be included with the help of field theories j^l- We note also that in our formalism it is 
also possible to study the entanglement of polymers in confined geometries. For example, 
values of E which are near to {E < a^^) correspond roughly speaking to the situation 
in which polymers fiuctuate in a quasi two-dimensional environment, in which the height in 
the t— direction is of the order of a few molecular sizes. 

Finally, an open problem, which has not be discussed here because we were mainly 
interested in the second moment, is the derivation of the full partition function Qx{E; ri, Tq) 
of Eq. (fT9|) . As anticipated in the Introduction, it is not an easy task to compute Gx{E; ri, Tq) 
because the repulsive potential of Eq. ^ is not central. We note however that the expression 
of Q\{E] ri, ro) coincides with the Green function of a spin 1/2 Aharonov-Bohm problem in 
the imaginary time formulation of quantum mechanics. This Green function has been already 
derived in [3| using sophisticated techniques developed in Refs. j^iSli which bypass all the 
difficulties of dealing with a non-central potential. Thus, in principle, the expression of the 
partition function ^a(-£'; ri, ro) is known. Unfortunately, some of the consistency conditions 
imposed on the parameters in Ref. [3| seem to be incompatible with the requirements of 
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our physical problem, as noted in Section UTTl For these reasons, the computation of the full 
partition function Q\{E; ri, Tq) is still a problem which needs further investigations. Luckily, 
the knowledge of the partition function is not necessary if one is interested to study the 



excluded volume interactions, which arise in the strong coupling 



imit. In fact, in this case 



it is possible to apply a powerful method due to Kleinert l^- This method turns 

the weak coupling expansion into a strong coupling expansion which is convergent for large 
values of vo and is able to accommodate also the anomalous dimensions of quantum field 
theories. The convergence of this strong coupling expansion is mostly very fast, so that only 
a few coefficients of the weak coupling expansion must be known, see Refs. for more 

details. These coefficients can be easily computed starting from the well known partition 
unction of the Aharonov-Bohm problem without the insertion of the 5— function potential 
3| and treating this potential as a small perturbation assuming that the value of Vq is 
small. The application of Kleinert's method in order to complete the brief analysis of the 
strong coupling limit made in this paper is work in progress. 
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